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Abstract We perform molecular dynamic (MD) simu- 
lations of frictional non-thermal particles driven by an 
externally applied shear stress. After the system jams 
following a transient flow, we probe its mechanical re- 
sponse in order to clarify whether the resulting solid is 
'fragile'. We find the system to respond elastically and 
isotropically to small perturbations of the shear stress, 
suggesting absence of fragility. These results are inter- 
preted in terms of the energy landscape of dissipativc 
systems. For the same values of the control parameters, 
we check the behaviour of the system during a stress 
cycle. Increasing the maximum stress value, a crossover 
from a visco-elastic to a plastic regime is observed. 

Keywords jamming • kinetic arrest • granular 
materials • rheology ■ phase diagrams 

1 Introduction 

Among the peculiarities of granular systems [JJ, there 
are jamming processes consisting in a sudden transi- 
tion from a flowing to a solid state [5J|3]. In particular, 
if submitted to an anisotropic stress, a granular mate- 
rial can jam after a transient flow. This transition is 
observed neither at low densities, where granular ma- 
terials behave like fluids, nor at very high densities, 
where they behave as disordered solids. The jammed 
solids obtained under these conditions, are so uncon- 
ventional that it has been claimed to represent an ab- 
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solutcly new state of matter, which have been termed 
'fragile matter' [3]. Fragility is based on the idea that 
during the flow particles rearrange until they find a con- 
figuration able to balance the external stress. It is spec- 
ulated that the force chains sustaining this configura- 
tion are strictly related to the flow-generating stress, so 
that they can support incremental stress oriented as the 
former (compatible stress). By contrast, the may not 
sustain a stress applied in a different direction (incom- 
patible stress), which may allow the system to restart 
flowing. This resembles the case of a pile of sands ob- 
tained in the presence gravity. The pile responds as a 
solids to stress directed as the weight, but it falls down 
under the action, for example, of a shear stress. 

A more precise statement of the idea of fragility has 
been provided in the limit of hard spheres [UG3E] • Con- 
sider, for example, a system jammed in the presence of 
a constant shear stress a zx : in this framework compat- 
ible stress variations, a' = aa zx with a > 1, lead to 
an elastic response. Conversely incompatible stresses, 
such as a' = aa zx + ba zy restore flow, also for infinitesi- 
mal values of b. In this sense, these systems are 'fragile' 
and strongly differ from any ordinary visco-elastic or 
elasto-plastic material. 

Previous results suggest that, in the presence of a 
shear stress, a fragile behavior can only be observed 
in frictional granular systems [ZUHMSJUnUIIlIISlIISlIISlIISJ 
[To] . Indeed, as reported in Rcfs. [TUB], flowing friction- 
less systems have never been observed to spontaneously 
select a microscopic state able to sustain the applied 
stress, while, conversely, in the presence of friction a 
much more complex phenomenology has been found [9j 
fTnl[TTl[T^[T8lpll[T5l[T6l . In particular, in recent works [H 
116) reporting molecular dynamics (MD) simulations of 
soft frictional grains at constant volume and applied 
shear stress, it was observed the phenomenology who 
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inspired the idea of 'fragile matter' [TJJ. Indeed, for 
some values of the control parameters, the system was 
found to select a configuration able to sustain the ap- 
plied stress after a small slip, or even after very large 
strains at constant rate [18]. 

Under these conditions, whether the mechanical re- 
sponse of the system is peculiar or rather similar to 
that of other amorphous and soft materials [19] is still 
an open question. In particular it is not clear to what 
extent fragility can capture the physics of these sys- 
tems. In more concrete terms, can small incompatible 
stresses restore the flow? 

In this paper, we investigate via MD simulations 
the mechanical response of frictional granular systems 
jammed in the presence of a shear stress. After shortly 
reviewing the numerical model and the overall phe- 
nomenology (Sec. [5]), we discuss the limit of validity 
of the concept of 'fragile matter' (Sec. [3]). We find that 
our system is not fragile but it responds as an elas- 
tic and isotropic solid to small incompatible stresses. 
Anisotropy and macroscopic rearrangements only emerge 
in the response to strong stress variations. In Sec. 2] we 
will focus on the stress-strain curves during a stress cy- 
cle, in order to check analogies with ordinary rhcologi- 
cal behaviours. Increasing the maximum stress reached 
during the cycle, we find a crossover from a visco-elastic 
to a plastic regime. 

2 Investigated System 

We perform MD simulations along the line of Refs [HI 
[TfilllSj . Monodisperse spherical grains of mass m and 
diameter D are enclosed in a box of dimension l x = 
l y = 16-D, and l z = 8D. Periodic boundary conditions 
are used along x and y, while the size of the verti- 
cal dimension is fixed and chosen to be comparable to 
that of recent experiments |20[|21| . The upper and lower 
boundary surfaces of the box are disordered collections 
of particles that move as a rigid object. The bottom 
plate has an infinite mass, and is therefore fixed, while 
the top one has a mass equal to the sum of the masses 
of its particles (roughly l x l y ). We impose a shear stress 
&zx — c to the system, applying a force to the top plate. 

Grains interact via the standard linear spring-dashpot 
model. Two particles i and j, in positions and r^, 
with linear velocities and Vj, and angular velocities 
uji and ojj, interact if in contact, i.e., if the quantity 
Sij = D — |ry| is positive. Sij is called the penetra- 
tion length, and ry = — Vj is the distance between 
particles i and j. The interaction force has a normal 
component F nij and a tangential one F tij : 

= -k n 8ijYiij - j n m e ffV nt: . (1) 
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Fig. 1 (Color online) Jamming phase diagram in the density- 
shear stress (<f) — <j) space at a fixed value of the friction, 
/i = 0.8. Three transition lines bound the different regions. 

F Uj =-k t u Uj . (2) 

Here k n and kt are elastic moduli, = r.y/|rjj|, 
v «, 3 - = [(vt - Vj) ■ nyjiiij, v t y = Vij - v n!j , m eff is the 
reduced mass, and 7„ accounts for dissipativc character 
of the normal component. Ut y , set to zero at the be- 
ginning of a contact, measures the shear displacement 
during the lifetime of a contact. Its time evolution is 
fixed by v tiJ and u>i and u>j, as described in Ref. [22) . 
The presence of tangential forces implies the presence 
of torques, = — l/2r.y x F ti . . We use the value of the 
parameters of [22]: k n = 2 10 5 , k t /k„ = 2/7, j n = 50. 
Length, masses and times are expressed in units of D, 
m and y/m/k n . We vary the volume fraction 0, which 
represents the volume occupied by the grains divided 
by the volume of the container, by changing the num- 
ber of particles. The initial state is prepared setting to 
zero the friction coefficient [25) , randomly placing small 
particles into the system, and then inflating them until 
the desired volume fraction is obtained; such a protocol 
is a short-cut of experimental procedures with which it 
is possible to generate very dense disordered states of 
frictional systems, such as oscillations of high frequency 
and small amplitude [2rJ] . 

We have recently investigated the jamming proper- 
ties of this system [16j as 0, <r, and /i are varied, and 
summarize our findings below. Confirming previous re- 
sults [7], in the frictionless case the system is either 
found in a flow or in a jammed phase, a single line 
marking the transition between these two phase in the 
<p-a (fx = 0) plane. By contrasts, at a finite value of the 
friction coefficient, more complex rheological regimes 
are found, corresponding to four distinct regions in the 
4>~a (/i > 0) plane, as illustrated in the schematic jam- 
ming phase diagram reported in Fig. [T] 

On increasing the volume fraction the system vis- 
its the following regimes: Flow, Flow & Jam, Slip & 
Jam, Jam. In each of these regime the system exhibits 
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Fig. 2 Typical behaviour of the velocity of the top plate as a 
function of time in the different regimes, from Ref. | 16| . Here 
a = 0.5, n = 0.8 and <j> = 0.578 (Flow), cj> = 0.596 (Flow & 
Jam), <j> = 0.629 (Slip & Jam), </> = 0.655 (Jam) 

a different response to the applied shear stress a, as 
illustrated in Fig. [5] where we show the time evolution 
of the velocity of the shearing top plate. The different 
regimes can be characterized as follows: 
Flow: the system flows with a steady velocity reached 
after a transient. 

Flow & Jam: the system reaches a steady velocity af- 
ter a transient. However, after flowing for a time, tj am , 
it suddenly jams. 

Slip 8z Jam: the system jams after a small inelastic 
displacement of the top plate. Steady flow is never ob- 
served. 

Jam: the system responds as an elastic solid. 

The value of the volume fraction (4>j 11 4>j 2 an d <pj 3 ) 
[27] marking the transition between the different regimes 
depend both on a and fi, and their identification [16| 
allow to describe the overall phenomenology in a three- 
dimensional jamming phase diagram in the 0, er, /i space. 
In the limit a — >• the transition from the 'Flow & 
Jam' to the 'Jam' regime appears to occur at a volume 
fraction identifiable with random-close packing volume 
fraction, if one neglects its protocol dependence [23] . 
We prefer not to associate the other transitions volume 
fraction to the random loose volume fraction, consid- 
ering that this has been introduced using a different 
protocol (pouring), and that this is also expected to be 
protocol dependent 

We have previously investigated the dynamical and 
geometrical properties in the different regimes |16j . For 
example, we find that in the 'Flow & Jam' region the av- 
erage value of tj am diverges on approaching the 'Flow' 
regime, thus allowing to identify the transition line <pj 1 
between 'Flow' and 'Flow & Jam' region. In addition, 
the tjam distribution is found to be extremely broad. 
Moreover, for each value of the friction, fj,, we find the 
mean contact number Z to be constant across the 'Flow 



Fig. 3 (Color online) a) Network of contact forces in a sys- 
tem jammed under the effect of a constant shear stress. Each 
segment marks the direction of the force acting between a 
pair of grains in contact. A color scale is used to represent 
the force intensities, b) The picture only shows the strongest 
contacts (10% of all contacts) of the network. 

& Jam' and the 'Slip & Jam' regimes. This is in agree- 
ment with Ref. [25] , even though we also found the con- 
stant value of Z to depend on the applied shear stress. 
We described in Ref. [TS] the rheology of the system and 
a possible mechanism for the transition from a flowing 
to a jammed state. 

In this paper, we fix the Coulomb friction coefficient 
to fi = 0.8, and consider values of <f> and a correspond- 
ing to the 'Flow & Jam' and to the 'Slip & Jam' regime, 
where the systems jams after flowing, and could there- 
fore be fragile. In the 'Jam' region fragility is not ex- 
pected as jamming is not preceded by a transient flow. 

3 Mechanical response of jammed states 

A minimal model for the structure of fragile systems 
consists in a scries of force chains directed along the 
strongest stress direction, supporting the applied stress 
and living in a sea of spectators [H[5j|6]. We show the 
network of all forces of a system jammed under the ac- 
tion of a shear stress in Fig. [5^. A percolating cluster of 
contacting particles furnishes the support for the net- 
work of force. Finite clusters, except few single parti- 
cles, called rattlers, are not allowed, since, due to repul- 
sive force acting between particles in contact, a cluster 
not hold by the confining plates breaks. The network 
may be also explored investigating how it changes as 
a function of a threshold on the inter-particle force 
one may introduce to define the bonds. For instance, in 
Fig. 13b we show the network of the strongest forces, ob- 
tained using a threshold which allow to pick-up roughly 
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Fig. 4 Response of a jammed system to a small perturbing 
shear stress. Panel (a) shows the response at a = 10 -2 , and 
8a = 10 -4 for different values of the volume fraction (from 
the inside, <f> = 0.655, 0.630, 0.617, 0.613 and 0.610). Panel 
(b) shows the response at <f> = 0.617, a = 10 -2 , for different 
values of the perturbing stress (from the inside, 5<j = 10 , 
5 10" 3 , 10" 2 , 2.5 10" 2 , 5 10" 2 , 7.5 10" 2 , 10" 1 ). 



10% of all forces (strongest ones). While the network of 
Fig. [3^. appears to be isotropic, that of Fig. [HJa appears 
highly anisotropic, with force chains lying in the di- 
rection of the strongest stress, and resembles that men- 
tioned by Cates et al. [HOE] ■ If the response to external 
perturbations is dominated by these strong forces, one 
may therefore expected a 'fragile' behavior. 

In order to check this point, we probed the elastic 
properties of a system jammed under the action of the 
existing shear stress, a zx = a, by superimposing a per- 
turbing shear stress. The non-zero components of this 
perturbing stress are 5a zx and <5cr 
that 8a 2 x + 5a 2 
therefore conveniently expressed in terms of 5a and of 
9 = axct&n(8<j Z y 1 5<J zx) ■ 

Figure 0] shows the displacement Sr = (5x, 5y) of 
the top plate position for different values of the vol- 
ume fraction at fixed a and 5a (left), and for different 
values of 5a at fixed a and <fi (right). Each curve is ob- 
tained applying a perturbing shear stress with (8 = 0), 
and then increasing 9 from to 2tt. While open path 
were expected for restored flow, we find that each curve 
describes a close path, which implies that the system 
responds elastically to the applied force. Moreover, this 
path is to a good approximation a circle (\5r\ ~ const), 
which implies that the elastic response is the same for 
all values of 9. An estimation of the degree of anisotropy 
in the response if obtained by measuring the parameter 



zx auu uu zy , we fix in such a way 
— 5a 2 . The perturbing shear stress is 



[5x 2 (9)+5y 2 (9)] 



5r 



1/2 -Tr 



where 5r = ([5x 2 (9) + 5y 2 (9)] 1/2 ) e . 

As illustrated in Fig. [5], the anisotropy of the system 
is small, being \£,(9)\ < 4% and does not reveal any 
particular pattern. 



(3) 




Fig. 5 (Color online) Anisotropy in the response to a small 
perturbing shear stress of a system jammed under the action 
of a large shear stress. Different curves refer to different values 
of the volume fraction. 
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Fig. 6 Using the same protocol of Fig. 0] response of a 
jammed system at <f> = 0.606, a = 10 -2 to a perturbation 
(5cr = 10 -4 . The open path signals the presence of restored 
flow. 



Being the response to small 5a elastic (the strain 
is proportional to the stress) and to a good approxi- 
mation isotropic (the strain does not depend on 9), as 
clarified by Fig. H] and Fig. [5l the system is character- 
ized by a well defined shear modulus G = lim^ _>o Safe, 
where e = 5r/L z is the shear strain induced by 5a. As a 
further characterization of the mechanical properties of 
the system, we have studied the Hessian along the line 
of Rcf . [55] , only finding zero eigenvalues (expected due 
to the presence of rattlers), and negative ones, while 
positive values are expected for fragile solids. 

This analysis clarify that frictional granular systems 
jammed at constant volume and applied shear stress 
are not fragile, as they respond elastically to small per- 
turbations, regardless to their orientation. The absence 
of fragility can be rationalized in terms of the prop- 
erties of the energy landscape of the system. Indeed, 
fragile jammed systems can be associated with saddle 
points, as their elastic energy may increase or decrease, 
depending on the direction of the perturbation, respec- 
tively leading to an elastic response or to an instability. 
Since dissipativc systems do not spontaneously arrest 
in an unstable point of their energy landscape, we ex- 
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Fig. 7 (Color online) Displacement of the top plate in a 
stress cycle. The stress is first increased to its maximum value 
cr m and then decreased to zero. Here a m = 5 10 _ 4 and /i = 
0.8. From right to left, <t> = 0.6475, 0.6477, 0.6480, 0.6482, 
0.6488. 




Fig. 8 (Color online) For a fixed value of the shear stress 
(<7Vn = 5 10 ), the residual strain e r decreases on increasing 
the volume fraction, and vanishes at a volume fraction (j>j 3 , 
which depends Oil (Tfn and /i. The straight line is a power law 
e r = a(<t> j3 - (j}) b , b ~ 1.2 and <j) j3 ~ 0.6495. 



pect them to arrest in a true energy minimum. Sys- 
tems that jam under the action of an applied stress 
are, therefore, not expected to be fragile. Of course, 
restored flow predicted by fragility, may appear in re- 
sponse to large stress variations, which are able to carry 
the system away from an energy minimum. Indeed, at 
higher da, we find that curves like those in Fig. [4] be- 
come more and more elliptic, signaling the emergence of 
anisotropy. At a threshold 8a c these curves turn in open 
paths, such as the spiral shown in Fig. [HI which evidence 
the presence of restored flow. The value of the thresh- 
old depends on the volume fraction, 8a c =5a c (4>) as well 
as on the applied shear stress. In the explored rage of 
volume fractions, we have never observed a fragile be- 
havior when the relative variation of the shear stress is 
Sa c /a < 1CT 3 . 

These results suggest that in the response to small 
perturbations the systems probes the roughly isotropic 
network of all forces of Fig. (3^,, rather than the strong 
network of Fig. \3]p. 



4 Stress cycle 

As discussed in the previous section, the mechanical 
response of jammed grains in the investigated range 
of control parameters resembles that of other elastic 
solids. However, granular systems are also expected to 
exhibit a plastic response, as other disordered systems. 
We have investigated this possibility probing the re- 
sponse of the system to stress-strain cycles. In [TB] we 
have used this approach to locate the line separating 
the 'Slip & Jam' and the 'Jam' region in the jamming 
phase diagram of Fig. [TJ Here we focus on the degree of 
plasticity emerging in the solid response for increasing 
value of the shear stress. 



We have measured stress-strain curves when the sys- 
tem undergoes a stress-cycle for values of the control 
parameters which span from the 'Flow & Jam' to the 
the 'Jam' region. Precisely, after preparing the system, 
we slowly increase the shear stress a up to a maximum 
value a m , and then decrease it to zero. After each stress 
cycle we measure the residual strain e = AL/L Z , where 
AL is the displacement of the top plate. 

At low volume fraction, in the 'Flow & Jam' regime, 
the initial state of the system is not jammed, and ac- 
cordingly we expect a finite residual strain at the end of 
the cycle, e r > 0. Conversely, at higher volume fraction 
the system has a solid response, we expect to be elastic 
(e r = 0) at small a m , and plastic (e r > 0) for a m over- 
coming the yield stress. Figure [7] shows the strain as a 
function of the shear stress for a small value of the maxi- 
mum stress, a m = 10~ 4 , and different volume fractions. 
Accordingly to the expectations, the residual strain be- 
comes smaller as the density increases, and appears to 
critically vanish at <pj 3 , which has been determined via 
a numerical fit, as shown in Fig. [8j On increasing a m , a 
finite residual strain e r > is also found at higher vol- 
ume fractions. This is shown in Fig.[S] where we observe 
that e r vanishes on increasing <f> at small a m , while it 
decreases and then bends at larger a m . 

The crossover in the behavior of e r can be related to 
a change from a visco-elastic to a plastic regime, and 
may be explained focusing on the increase of the num- 
ber of contacts that break as the strain increases. At 
small stress, the strain of the system is small, contacts 
do not break and the system responds elastically. At 
higher stress, the strain of the system is large, and con- 
tacts break. This is the microscopic origin of the plastic 
response, as the system looses memory of the tangential 
force of the contacts reaching the Coulomb threshold. 
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Fig. 9 (Color online) Dependence of the residual strain e r on 
the volume fraction <f>, for different values of cr m , at fj, = 0.8. 



5 Conclusion 

We have investigated the mechanical response to small 
perturbations of systems jammed under the action of 
an applied shear stress, and found this to be elastic and 
isotropic, and therefore not consistent with the expecta- 
tion of a 'fragile' state. This result is explained consid- 
ering that the response to small perturbations see the 
cooperation of all contact forces, which are isotropically 
distributed in the systems. Conversely, the response to 
large perturbation is due to the strongest forces, which 
are not isotropic. Accordingly, we have found that large 
perturbations lead either to an anisotropic response, 
or to unjamming. The investigation of the response of 
the system to stress-cycle, also allowed to observe the 
transition from an elastic to a plastic response as the 
maximum applied stress increases. 
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